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Abstract 



A generalized approach to the quantization of a large class of maps on a 
^ ■ torus, i.e. quantization via the von Neumann Equation, is described and a 

■ number of issues related to the quantization of model systems are discussed. 
The approach yields well behaved mixed quantum states for tori for which the 

m 

On . corresponding Schrodinger equation has no solutions, as well as an extended 

J>^. spectrum for tori where the Schrodinger equation can be solved. Quantum- 

Q . classical correspondence is demonstrated for the class of mappings considered, 

■ with the Wigner-Weyl density p{p, q, t) going to the correct classical limit. An 
application to the cat map yields, in a direct manner, nonchaotic quantum 

Ix! 

'Vl . dynamics, plus the exact chaotic classical propagator in the correspondence 

^. 

limit. 



I. INTRODUCTION 

Quantization of systems with classical analogs proceed by a well known procedure. 
Specifically, one replaces coordinates or momenta by operators in the Hamiltonian, con- 
structs the Schrodinger equation and imposes boundary conditions which arise directly from 
the classical physics. In this paper we show for the case of maps on a torus, model systems in 
nonlinear dynamics, that modifications to this procedure are essential to allow the quantum 
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dynamics to unambiguously approach the correct classical limit. 

Area preserving mappings on the torus are model systems which have long been used 
to illustrate the essence of regular or chaotic classical dynamics because they are of low 
dimensionality and can be readily propagated. These maps are essentially systems which 
are confined to finite coordinate and momentum space regions 0<g<a, 0<]9<6, with 
rules of dynamic evolution which confine the system to this region through a dependence 
on q mod a, p mod h. In addition, time evolution is in discrete steps. Much of the cur- 
rent interest in quantizing these mappings stems from the expectation that they will 
prove useful in understanding the role of classical chaos in quantum dynamics and in under- 
standing quantum-classical correspondence. Indeed, little is known about quantum-classical 
correspondence for systems whose classical analog is chaotic [0. While the correspondence 
principle requires that quantum mechanics should produce the classical laws in the limit 
that Planck's constant approaches zero §], the characteristics of the limit are such as to 
make this verification extremely difficult ||^. 

Central to the issues of classical-quantum correspondence is the specific prescription 
used to quantize a given classical system. In this paper we show that a new extended 
quantization procedure allows an alternative method of quantizing mappings on a torus 
which yields solutions having a number of desirable features which do not arise in the 
traditional quantization via the Schrodinger equation. 

Motivation for reconsidering the quantization procedure for maps on a torus stems from 
an examination of the results of previous efforts in this area. Hannay and Berry p[, and 
Balazs and Voros 0, chose to directly quantize the Schrodinger propagator while Ford 
and coworkers [|l| quantize certain maps by introducing a kicked oscillator Hamiltonian 
{—oo < q,p < oo) 

oo 

i/ = pV2/x + egV2 S{s-t/T), (1) 

s=— oo 

The dynamics generated by this Hamiltonian stroboscopically produces the classical dynam- 
ics of the map. They then construct the propagator for this system. The restriction of the 
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classical dynamics to a torus 0<g<a, 0<j9<6is introduced by imposing periodic bound- 
ary conditions on the wavefunctions in both the position and momentum representations. 
The results of this procedure have been the subject of considerable controversy. Specifically, 
Ford 0] showed that although classical cat map dynamics is algorithmically complex, the 
quantum cat does not show this type of random behavior in the classical limit. Berry and 
coworkers, however, argue 0] that such measures are less useful indicators of correspon- 
dence failure than believed by Ford and coworkers. Further, Graffi et al |^ have recently 
shown that the quantum cat is mixing in a specific limit, N ^ oo along the subsequence of 

prime, where h = 1/N. This is in contrast with Ford's demonstration of quasiperiodic 
dynamics when taking a somewhat different approach to the classical limit. 

Clearly the traditional quantization procedure leads to a quantum cat map whose be- 
havior as /i — > is, at best, subtle. In addition the traditional approach leads to two other 
features which we regard as undesirable: (a) only a restricted class of classical tori, those 
with ab = hN, N integer, can be quantized ||T0[. Indeed, some mappings can only be quan- 
tized for even and (b) the resultant wavefunctions are highly singular, being a set 
of delta functions. 

In this paper we show that introducing a translationally invariant Hamiltonian, quan- 
tizing the associated von Neumann (quantum Liouville) equation dp/dt = [H,p]/ih for the 
density matrix p, and obtaining solutions satisfying Bloch boundary conditions, allows quan- 
tization of tori for all ab, is physically complete, leads to smooth well behaved solutions and 
is measure preserving. Further, we show that the correct classical propagator and the exact 
classical dynamics emerge directly from the quantum propagator in the classical limit, with 
the limit taken as the complete sequence N oo. Hence the classical limit of our quan- 
tum propagator for the cat map generates algorithmically complex dynamics, the system is 
mixing, etc. 

Note that our approach yields density matrix solutions for cases (e.g., ab ^ hN) for which 
there are no pure states. The fact that density matrix solutions exist when wavefunctions 
do not is of considerable interest and is emphasized below. 
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This paper is organized as follows: Section II summarizes aspects of the traditional 
quantization procedure and treats the case of e = in Eq. (1) explicitly. The extended 
quantization procedure is then introduced in Section III where it is also applied to the case 
of Eq. (1) with e = 0. Examining this case allows us to emphasize the different spectra and 
solutions obtained by this approach and to show that both eigenfunctions and eigenvalues of 
the von Neumann equation directly approach the classical limit. In Section IV we generalize 
the discussion to the case of e 7^ and use the von Neumann quantization scheme to obtain 
the quantum propagator. Classical and quantum map dynamics are expressed in a form 
suitable for the study of correspondence. Section V then explores the quantum dynamics 
in the correspondence limit where we show that all maps treated reduce to the correct 
classical limit. In the application to the quantum cat we show that its quantum dynamics 
is not ergodic, but that the chaotic classical cat is nonetheless recovered in the h ^ 
limit. Throughout this treatment we make extensive use of the Wigner-Weyl representation 



11,12], a p,q dependent representation of quantum mechanics whose form encourages a 



direct comparison with classical dynamics. In Section VI we address the question of why 
there are no pure states for the quantum cat, showing that the system best resembles a 
discretized Langevin process. Finally, Section VII contains a summary and remarks on 
future work. 



II. TRADITIONAL QUANTIZATION 

Constructing the appropriate quantum generalization of a given classical system requires 
a correspondence rule for replacing classical observables by quantum operators, a choice of 
a dynamical state equation and a statement of the associated boundary conditions. The ap- 
propriate dynamical state equations are normally assumed to be the Schrodinger equation 
for pure state dynamics and the von Neumann equation (or Quantum Liouville Equation) 
for mixed state dynamics and, customarily, the choice of boundary condition is dictated 
by features of the physical system. Quantum mappings on the torus have no known phys- 
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ical realization, and so some freedom would appear to exist for the choice of boundary 
conditions. Below we discuss the limitations of quantization via the Schrodinger equation, 
both for (previously adopted) periodic boundary conditions as well as for Bloch boundary 
conditions. We later show that a physically complete theory of quantum mappings on a 
torus does result from a von Neumann based quantization procedure in concert with Bloch 
type boundary conditions. Of general interest is the different way in which the boundary 
conditions, implemented in these two procedures and both apparently consistent with the 
classical map, affect the quantum picture. 



A. Boundary Conditions 

Traditional quantization procedures define the dynamics on the full phase space, 
— oo < p,q < oo, but impose periodic boundary conditions on the wavefunction in both the 
coordinate [V'(q')] and momentum representations [^/'(p)], i.e. 

^Piq) = ^|Jiq + a)=T^ia)^^Jiq) (2) 

i^{p)=i^{p + b)=Tp{bmp) (3) 

where Eqs. @ and (|^) implicitly define the coordinate and momentum translation operators 

Tq{a) = exp{ipa/h)] Tp{h) = exp{—iqb/h). (4) 

As a consequence of Eqs. @ and @, dynamics within each {q,p) unit cell is expected to 
reflect the character of torus dynamics [i.e. dynamics with variables (q mod a) and (p mod 
b)]. From Eq. (^) it follows that {1 — e~^^'^^^}'ip(q) = which implies that ip{q) = unless 
q = nh/b where n E Zi. Similarly, from Eq. (|^) it follows that {1 — e^'^^^^}'(jj{p) = which 
implies iplp) = unless p = mh/a where m G Z. Imposing both conditions simultaneously 
gives {1 — e*°*''^}'?/'(g) = and {1 — e~*''^/^}'?/^(p) = so that we must have ab = hN where 
A G Z. Fixing A one then obtains the class of allowed wavefunctions 
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^(g) = v^ E (5) 

j=-oo 

with ipj+N = V'i- This condition imphes that 

N 

^,=Y.c^,e'^''"\ (6) 
1=1 

where the ai coefficients are arbitrary. In the momentum representation these states take 
the form 

oo TV 

m = ^bY.T. ^kSip - (j + k/N)b) (7) 

j=—oo k=l 

where ipk = VNak- Thus, in general (independent of the Hamihonian H) periodic bound- 
ary conditions in p, q lead to a discrete spectrum with non-normalizable wavef unctions, 
comprised of sets of delta functions. The traditional interpretation of {ipl'^ as a probability 
is inapplicable. Note that it is impossible to quantize any mapping on a torus by this route 
unless ab = hN, a consequence of the Fourier transform relationship between p and q and of 
the requirement to simultaneously fit waves in both the q and p directions. With ab ^ hN 
neither the Hamiltonian, nor any other operator, has eigenvalues. The quantization proce- 
dure is thus incomplete insofar as only a small number of classical tori can be quantized 
||T0| . Finally, note that these conclusions are independent of the system and so apply to any 
mapping confined to a torus. 



The Wigner-Weyl picture |]TT],|T2| provides a p,q dependent quantum representation in 
the density matrix formulation which allows a direct comparison with classical mechanics in 
phase space. One can readily demonstrate that in this picture the dynamics lies on points 
in phase space which are rational multiples of a and b. Ford shows |l| that this behavior 
survives in the classical limit and is responsible for non-algorithmically complex dynamics 
for the classical limit of the quantum cat. It is obviously also the case for any map quantized 
by this procedure. 

The situation is similar if Bloch type boundary conditions are adopted. That is, Bloch 
boundary conditions on the wavefunction 
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^{p + b) = e-''^'>^/^^{p) (8) 
with constant qo, po again yield states 

oo 

H<l) = E -<1o- Ja/N) (9) 

j=-oo 

only when ab = hN. Here ipj = e^po^^i/^^ Y^^=i a^e^'^*-''/^ where the ai coefficients are arbi- 
trary. In the momentum representation 

oo TV 

^(p) = v^e-*^"?°/" E E ^kKv -Po- U + k/N)b) (10) 

j=— oo k=l 

where ipk = VNa^- Bloch boundary conditions will, however, prove useful in the extended 
quantization scheme discussed below. 



B. The case of e = 



Consider now, for illustration purposes, the case of e = 0. Here Eq. (1) gives = ab/h 
independent Hamiltonian eigenfunctions ||13|| " 



ij.ip) = ^b/N y: Sip - u + i/m (11) 

j=-oo 

and associated discrete spectrum Ei = {Ib/Ny /2fi, I = 1,...,N. We also consider, for 
comparison with results below, the solution from the perspective of the von Neumann (or 
quantum Liouville) equation: 

dp/dt = -i[H, p\/n = -iLgp. (12) 

Here p is the density operator and Eq. ([T2|) defines the quantum Liouville operator Lq. The 
eigenvalues A^^; and eigenfunctions p^^i of the Lq satisfy: 

LgPk,i = j:[H, pk,i] = \k,iPk,i (13) 

Both may be obtained directly from the Hamiltonian solutions and the latter conveniently 
expressed in the Wigner-Weyl representation. Specifically, {k,l = 1, ...,N): 



Xk,i = [Ek - Ei)/h = [k' - l')b'/{2f,N') (14) 

and pk,i = {ipk >< 'ipil so that the Wigner-Weyl representation of pk,i, denoted Pk,i{p,(l), is 
given by |11],|12 

/oo 
dq'<q~ q'/2\^k >< + q' /2 > e^^« 
-oo 

/oo 
dp' <p + p72|^fc ><iji\p- p'/2 > e'^P /^ (15) 
-oo 

Using Eq. (pUf) we have 

p^i{p, q) = 2a-^ e2"*("-'")'?^/"e2"'('=-')«/"5(2p -{n + m+ -—)b). (16) 

n,m=— oo ^ 

A simple change in the integers of summation in this expression for pk,i{p, q) and use of the 
identity m 



^ 6^-^^-^= E 5{x-3) (17) 



j=— oo j=— oo 

gives 



oo 



Pfc,KP,?) = (2iV)"'e2-('="'^<?/'^ E (-l)™"%-ma/2iV)5(p-(n + ^)6/2). (18) 

n,m=— oo 

Thus, the eigenfunctions of Lq in this representation are delta functions localized on a 
discrete number of points of the torus which are rational multiples of a and 6, i.e., a delta 
function "brush" in phase space. Note that Eq. ([17|) is repeatedly used throughout this 
paper although it is not always cited. 

The results in this subsection are used later below for comparison with the extended 
quantization results. 



III. EXTENDED QUANTIZATION 

An examination of a different approach to map quantization procedure is hence well mo- 
tivated. The extended quantization procedure which we advocate entails three new features: 
(a) the introduction of "mod operators" which produce a properly symmetrized Hamilto- 
nian, (b) solving the von Neumann [Eq. ([12|) ], rather than Schrodinger, equation and (c) use 
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of Bloch boundary conditions on the density matrix. The latter is motivated immediately 
below. [Throughout this paper, as above, we denote quantum density operators by p and the 
Wigner-Weyl representation by p{p, q). In addition, analog quantities in classical mechanics, 
such as the phase space density, will carry a "c" superscript.] 

A. Boundary conditions 

Consider then the von Neumann equation with either of the two choices of boundary 
condition i.e. periodic and Bloch type. From the perspective of the density matrix, periodic 
boundary conditions on the wavefunction imply Tp = p, pT = p for T = Tq{a) and T = Tp{b). 
In the Wigner-Weyl representation pJ| 

/CO 
dq' <q-q'/2\p\q + q'/2> e'P''^^ (19) 
-oo 

these conditions take the form 

e'P^/^pip, q ± a/2) = e-'P'^/^pip, q ± a/2) = p(p, q) 

e^'^'/^pip ± b/2, q) = e-^'^'/^pip ± 6/2, q) = p{p, q). (20) 

The following requirements on the density result from Eq. (^Dj): (a) (1— e*^''*/^)p(p, q) = Q 
which implies q = hl/2b where / G Z, (b) (1 — e*^P'^/^)p(p, g) = which implies p = hk/2a 
where G Z, and (c) (1 — e*"^/^)p(p, g) = which implies that ab = hN where G Z. 
The resulting states are mixtures of the delta type wavefunctions of Eq. (^) localized on the 
rational points of the torus and they exist only when ab = hN. Hence this combination of 
dynamic law and boundary condition regenerate the same difficulty as that associated with 
the Schrodinger approach and is, once again, independent of the system. This is the case 
since Eq. ( pOf ) is an exact statement of the requirement of periodicity on the wavefunction, 
transferred to the density matrix. It results in precisely the same restriction as it did in its 
application in the wavefunction picture 

A similar situation does not, however, arise with the application of Bloch boundary 
conditions to the density matrix p. Specifically, here Eq. (^) implies that 
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< g + a\pjk\q + a >=< q + a\ipj >< ipk\q + a >=< q\ipj >< ipk\q >=< q\pjk\q > (21) 

with a similar result in the momentum space picture. Thus, the phase exp{ipQa/h) asso- 
ciated with the wavefunction drops out, resulting in a quite different boundary condition 
in the density matrix picture than in the wavefunction picture. Specifically, the conditions 
T^^pT = p imply Bloch boundary conditions on the wavefunction ip ii p = but den- 

sities p can exist which satisfy T^^pT = p and which do not take the form of a weighted sum 
over pure states satisfying Bloch boundary conditions on the wavefunction. These boundary 
conditions will prove extremely useful in conjunction with the introduction of mod operators, 
as described below. 



B. Mod Operators 

We consider solutions to the von Neumann equation with Bloch type boundary condi- 
tions. To do so we first introduce the operators (pmod b) and (gmod a). These are obtained 
as a direct extension of the classical Fourier expansions for {q mod a) and {p mod b), i.e. as 
extensions of : 

1 °° i 
(gmod a)/a = Qmod 1 = - + V fo,m 

2 mt^oo 2vrm 

1 °° i 

(pmod = Pmod 1 = - + J2 7^fn,o, (22) 

where fn,m{P,Q) = exp[27ri(?2P + mQ)] and where we have introduced the convenient di- 
mensionless variables Q = q/a, P = p/b. Specifically, the quantum operators, analogs of Eq. 
(0), are 

^ 1 °° i ^ 

(pmod 6)/6 = Pmod 1 = - + ^ 7— /„,o, (23) 

with eigenvalues Pmod 1 and conjugate variable P!7[] ° 

^ ^ I °° i ^ 

(gmod a)/a = Qmod 1 = - + ^ /o,m (24) 

^ m=-oo ^T^fn 
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with eigenvalues Qmod 1. Here the unitary operators fn,m are obtained by Weyl quantizing 
H fn,m{P, Q) = exp[27rz (nP + mQ)], i.e. : 



fn,m = I dp I rfge2"("^+™^) I dve'P''% + v/2){q-v/2\ 

/oo 
dq e^"™^|g - naa/2){q + naa/2\ 
-oo 

/oo 
dp e2™-P|p + mab/2) {p - mab/2\ (25) 
-oo 



where a = h/ab is a dimensionless form of Planck's constant. Although Eq. ([25| ) provides 
forms most useful for computations, applying the coordinate translation operator inside the 
integral allows us to rewrite Eq. (^) in a far more attractive form: 

= exp{27ii{nP + mQ)}. (26) 

where (Q,P) = {q/ci,p/b) are the scaled coordinate and momentum operators. The fn,m 
operators whose properties are discussed in Appendix A, will allow us to write the 
quantum operator analog of the Fourier transforms necessary to treat periodic systems on 
tori. 

The introduction of these operators allows us to revise the Hamiltonian in Eq. (1) to 
read: 

oo 

H = {pi^dbf/2fx + e{qir^daf/2 ^ 5{s-t/T). (27) 

s=— oo 

This new Hamiltonian satisfies [H, Tp{b)] = [H, Tg{a)] = 0, i.e. H reflects the desired cell-like 
character of the full phase space. This character is also demanded of the eigenfunctions pij 
of the von Neumann equation, i.e. by requiring the generalized Bloch boundary conditions 
on p, 

p = T~\b)pT,{b) 

p = T-\a)pTg{a). (28) 
These conditions have the form, in the Wigner-Weyl representation, of: 
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pij), q + a) = p{p, q) 

p{p + b,q) = p{p,q). (29) 

Hence these conditions correspond precisely to the boundary conditions one would impose 
on a classical phase space density if the density (rather than the wavefunction) were periodic 
in p, q. 

States satisfying Eqs. (p9|) must, by Fourier's Theorem, be of the form 

OO oo 

p{p,q) = ^ Y1 Pn,mfn,m = ^ Pn,mexp{2ni{nP + mQ)} (30) 

n,m=— oo n,m=— oo 

in the Wigner-Weyl representation and so distributions satisfying Bloch boundary conditions 
are expected to satisfy the analogous quantum Fourier expansion: 

oo oo 

p=— Pn,mfn,m = —r Y Pn,m exp{27ri (nP + mQ) } . (31) 

n,m=—oo n,m=—oo 

We utilize this approach below to examine the case of e = 0. Specifically, we solve the von 
Neumann equation [Eq. ([1^)] using the Hamiltonian of Eq. (pTj) coupled with the Bloch 
boundary conditions which require Eq. (|31|). The general case e 7^ is considered in Section 
III. 

C. The case of e = 

Consider the extended quantization procedure applied to the case of e = 0. The relevant 
Hamiltonian is then ||13|| ^ 



H = (pmod 6)V2/i = ^^(Pmod if /2p (32) 
and the associated von Neumann equation 

Wp = ^[(^mSd l)\ p] = idp/dt. (33) 
Eigenf unctions of the Liouville Operator satisfy: 

Wp = ^p- (34) 
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Inserting Eqs. (pTD and (|3^) into Eq. (|3^) we obtain 

Jj2 oo ^ ^ oo ^ 

Pn,m [(^'mod 1)^ = A ^ 

Pn,mfn,m- (^^) 



?i,m=— oo ri,m=— oo 



Using Eq. (ESf) it can be readily shown that, for arbitrary Pq (i.e. arbitrary po/b), 



oo 

E e-'™''1(Pm^dl)^/„,^J 

?i=— oo 

oo 

= {[(Po + ma/2)modl]2-[(Po-W2)modl]2} ^ e-2™^V„,^. (36) 

n=— oo 

Hence we can choose pn,m = ^m,jG~'^''^"'^° ■ The eigenfunctions and eigenvalues of Lg with 
Bloch boundary conditions are then given by: 

oo 

p^.p^ = (a6)-i E e2"*["(P-P«)/''+^'^/'^l (37) 

n=— oo 

or, alternatively, using the momentum representation [Eq. (^Sj)] of and Eq. (^) 

oo 

Pi.po = E bo + n6 + vrj/i/ a) {po + nb - njh/ a\ (38) 

n=— oo 

with eigenvalues 



where j is an integer, < po < b, and the arguments of the bras and kets in Eq. (|38D 
correspond to eigenvalues of p. The Wigner-Weyl representation p^^pg {p, q) of these solutions 
makes evident their character: 

oo 

P,AvA) = a~^e^^'^'^''^ E 5{p-po-nb). (40) 

ri=— oo 

Stationary states {\j,pQ = 0) are of two types, arising from either j = or po = 0. The 
stationary states with j = are uniform in q and belong to the point spectrum, whereas 
the p, q integral over the stationary states with j 7^ and Po = are zero, characteristic of 
elements of the continuous spectrum. The overall spectrum of Lg is continuous, due to the 
continuous po label. 
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Several aspects of these solutions are worth emphasizing. First, Eqs. (^) and (|39D 
provide solutions for all values of a and h. Second, the eigenfunctions in Eq. ( |38D are, in 
general, mixed states, since wavefunctions exist only for ah = hN whereas Eq. ( pHD applies 
to all ab. If the condition ab = hN is satisfied then it is possible to recombine some of the 
degenerate Liouville eigenstates to produce pure states. Even then, however, the pure states 
are only a small fraction of all of the possible solutions. 

To permit a comparison of Eqs. (0), (^1) and (^) to the classical limit, consider 
solutions to the classical Liouville equation for motion on a torus with the appropriate 
classical boundary conditions [which are the same as Eq. (^Of)]. That is we consider the 
classical dynamics via the Liouville equation: 

5p(^) /dt = {H, p(^) } = -iL,p^'^ (41) 

Here {, } denotes the Poisson bracket, p^'^^ is the classical phase space density and this 
equation defines the classical Liouville operator L^. Eigenfunctions and eigenvalues of L^. 
are given by the solutions to the classical Liouville eigenvalue problem 

L^p(^) = -i{pmod b/p)dp^''>/dq = A("V^'^- (42) 

Solutions to Eq. (|4^ ) are expected to satisfy Eq. (pQ]). Hence p^^^ must take the form of Eq. 
(0) and so (with classical Fourier coefficients being denoted p^'^^^) 

n,m=-oo ' 

Using Eq. (pSf ) and taking matrix elements in the Fourier basis 



^ paJ^--- \a,h^k-n ^^^^ 

We assume a solution of the form p^^^^ = ^2^6"*"^"" and substitute this expression into Eq. 
(0). We obtain 

[AW - ^]l]„e-"^- = -,I^^]„e-^"^- f ^ (45) 



fc#0 
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or 



Using the identity 



we have 



nmb 2mb ^ sin(A;6'^ 



IJ,a 



k=l 



k 



(46) 



E 

k=l 



sm{k9) _n-9 
k ~ 2 



(47) 



On 



mb 



(48) 



Defining 0m/2vr = po/^, Po ^ [0, &), and choosing = Smj so that p^*^)^ 
gives resultant classical eigensolutions 



X -2mnp(i/b 



1 27rjjg/a 



E 



^27rm(p-po)/6 



n=— oo 

oo 



^-1^2nijq/a 



E 5{p~po-nb), 



(49) 



with associated eigenvalues 



(50) 



with < po < ^- Here the second identity results from an application of Eq. ([171). 

A comparison of the classical [Eqs. (^) and (0)] and quantum [Eqs. (0) and (p9D] 
results shows that the quantum and classical eigenf unctions are identical and the quantum 
eigenvalues go to the proper classical limit as h ^ 0. This is not the case with the results 
of the traditional quantization procedure, as discussed by Ford 



D. Remarks on Extended Quantization 

The procedure above involves two modifications to the standard quantization approach: 
(a) the introduction of modular variables and operators into both the classical and quantum 
problem and (b) the use of the von Neumann equation in quantum mechanics. The need 
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for the latter stems from the fact that the solutions which we find are, in fact, not pure 
states and are not comprised of averages over pure states. That is, we find that solutions 
to the von Neumann equation exist where Hamiltonian eigenf unctions do not. The need 
for modular variables and operators is also worth emphasizing. Specifically, note that the 
solutions Eq. ( P^ and Eq. (^D|) are not eigenfunctions of their respective quantum and clas- 
sical Liouville operators if these operators are written in terms of nonmodular p, q variables. 
Indeed eigenfunctions of the Liouville operators with Cartesian p, q variables and periodic 
boundary conditions do not exist for mappings on the torus. This is easily seen as follows. 
Consider the Liouville problem in p, q. Since the Hamiltonian in Eq. (1) with e = is 
quadratic it gives rise to a quantum Liouville operator which, in the Wigner-Weyl represen- 
tation, is identical to the classical Liouville operator . Thus, to obtain both the quantum 
and classical Liouville eigenfunctions for e = we would need only consider the classical 
Liouville eigenequation Lcp'^^^ = X^p^'^^ where Lc = —i{p/jj,)d/dq. Taking matrix elements in 
the Fourier basis would again yield Eq. which as we have seen has the solutions given 
by Eq. (^). But consideration of these solutions shows that they are eigensolutions of 
Lc = —i{pmod b/fi)d/dq, and not of Lc = —i{p/p)d/dq. Thus the procedure, starting with 
Cartesian p, q variables generates an inconsistency whose origin is readily apparent. That is, 
although the Hamiltonian is of the form given by Eq. (1), the Hilbert space chosen to solve 
the problem is that spanned by the Fourier basis functions. But the Hamiltonian in Eq. 
(1) and its associated Liouville operator cannot be decomposed on this Fourier basis, i.e., 
they do not satisfy the translational invariances of the Fourier basis. Hence, mathematical 
inconsistency and an incorrect form for the quantum propagator result. 

The advantages of the von Neumann approach over the conventional Hamiltonian quan- 
tization are now evident. Equations ( pSf ) and ( p9D provide a solution for any value of ab, are 
integrable over any of the periodic unit cells and are well behaved. Further, in the classical 
limit {h — » 0), Aj^pQ A^-^p^ [compare Eqs. (pQ]) and (|50|)] so that the classical limit is 
properly reached by both eigenfunctions and eigenvalues. 
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IV. CLASSICAL AND QUANTUM MAPPING DYNAMICS (e / 0) 



Here we focus on the more general case of e 7^ 0, obtain the propagator and show that it 
readily yields the correct classical limit. 



A. Classical Mechanics 

Consider then the Hamiltonian [Eq. (1)] in reduced modular variables [r] = Tb/fia e Z 
and ^ = -eTa/b G Z): 

ah 



^-^ s=-oo 

The classical dynamics of the map is defined by the classical propagator 

»(M+i)r 



(51) 



Ae(T) = rexp{-z / dtL,{t)} 

JMT 



(52) 



where L^if) is the classical Liouville operator associated with Eq. ( pTD and where the integral 
is over a time interval beginning just after the M*^ kick to just after the (M + l)*** kick. 
Then the dynamics is given by 



Ae(r) 



^ gmod 1 ^ 






^ Pmod 1 ^ 







mod 1 



(53) 



where 



(0 \ 


= 


Qn 


mod 1 = 


^ 1 7] ^ 


Qn 












[Puj 






[Pn) 



mod 1. 



(54) 



These linear mappings are such that I = ^Q^ + r]^QP - r]P^ is a constant of the motion 
when |2 + r]^\ < 2. However, when |2 + rj^l > 2, / is a hyperbola which is wrapped densely 
around the torus by the "mod 1" operation and so effectively destroyed. We examine the 
quantum dynamics of the map in both these regimes and obtain results for all rj, ^, including 
the Arnold cat map (77 = ^ = 1), and systems like the quasiperiodic discretized particle 
in a box (^ = 0). 
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To obtain the time evolution of the phase space density which at time zero is p^'^\P, Q, 0), 
we Fourier expand in fnm- At time zero, and later time MT, we have: 



-| oo 

p('=)(P,g,0) = - E Pit(0)/n,^ (55) 



n,m=— oo 



-| oo 

,(^)(P,Q,Mr) = - P^nUMT)fn,m = AciTfp^^\P,Q,0). (56) 



n,m,=— oo 



A straightforward analysis of the effect of A^T) on p^^\P, Q, 0) gives 

A,(r)p('=)(P, Q,0) = ^Y: Pt'inm)fn,m (57) 

n,m 

where (f)'{nm) is the transpose of the matrix in Eq. (|54D multiplying the column vector 
with elements n, m. Clearly, evolution of the map, from this perspective, corresponds to the 
interchange of coefficients in the Fourier basis fn,m- We define the classical propagator Gc 
of the Fourier coefficients, for particular reference with the quantum result obtained later 
below, as 



pi%iT)= Y: a(n,m;A;,/)pg(0). (58) 

fc,/=— oo 

Comparison of Eq. (^) for M = 1 and Eq. (0) reveals that 

Gc{n, m; k, I) = 6i^k,i),(i>'(nm)- (59) 
Generalization to the step from time MT to time (M + 1)T gives 

p(t((M + 1)T) = Y a(n, m; k, l)p'^i{MT)) (60) 

k,l 

SO that 

p('=)(P, Q, (M + l)r) = 4 EE Gc(n, m; fc, /)pg(MT)]A,^. (61) 

'n.,m k,l 

Equation (pl| ) provides a general expression for the time evolution of the classical phase 
space density in terms of the classical propagator Gc matrix elements as coefficients in a 
Fourier expansion. 



B. Quantum Mechanics 



Consider now the quantum dynamics with the Hamiltonian, the analog of Eq. (^) with 
the appropriate introduction of mod operators, 

ah °° 

H{t) = Tp^MP^od ir - e(gmod 1)2 Y: ^(^ - t/T)}- (62) 

-^-^ s=-oo 

The associated time evolution operator is 

U{T) = Te-'^MT^^^ "^'^^^l^ (63) 



p(MT) = -i- PnAMT)fn,m (64) 

ra,m=— oo 



where the integration is performed from just after the M*'^ kick to just after the (M + 1)** 
kick. The resultant U{T) is M independent. 

The goal of our calculation below is to obtain the quantum analog of Gc- To do so is 
complicated in execution but simple in its essence. Specifically, we introduce the Fourier 
operator basis [Eq. (126|) ] for the quantum distribution p at time MT as 

oo 

ab 

and endeavor to express the time evolution of the Fourier coefficients pn,m as the quantum 
analog of Eq. (|60|), i.e. as 

oo 

Pn,UiM + l)T)= G{n,m-k,l)pk,i{MT). (65) 

fc,i=— oo 

Doing so will allow an analysis of the quantum dynamics in the Fourier basis as well as a 
direct comparison with the classical result. Indeed this comparison with the classical limit 
benefits considerably from use of the Wigner-Weyl representation. In this representation 
Eq. ( p^ assumes the form 

1 oo 

piP,Q,MT) = - Y PnAMT)fn,m (66) 

„,m=-oo 

SO that, given Eq. (|65|), 

p(P, Q, (M + 1)T) = 4 EE G{n, m- k, l)puAMT)]fn,^. (67) 

n,m k,l 
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A comparison of Eq. (|67|) with Eq. ( |6lD shows that classical-quantum correspondence 
requires 
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YimpuAMT) = p^,{MT) (68) 

a— >U 

and 

\im.G{n,m]k,l) = Gc{n,m]k,l). (69) 

a— >0 



Equation (|g) holds if Eq. (^) holds and if Pk,i{^) = Pkli^)- However, the latter equality 
is assured, as seen by comparing Eq. (^) with Eq. (|55D , at time zero. Thus to prove 
the correspondence limit for the fundamental dynamical entity p{P, Q, MT) requires that 
we show Eq. (|69|). This is done in Section V after we obtain a general expression for the 
quantum G{n, m; k, I). Readers uninterested in the details of this derivation should proceed 
toEq. (II). 

1. Quantum Propagator 

The time evolved density operator satisfies 

p((M + 1)T) = tj{T)p{MT)U\T). (70) 

We transform to the Wigner-Weyl representation and focus on the initial step from time 
to T; the same algebra pertains to the general step from time MT to (M + 1)T. Equation 
([70|) in the Wigner-Weyl representation is 

p{P,Q,T) = - ^ p„,„ / dve'^-''''/"{q-v/2\KUo{T)fn,mU^o{T)K^\q + v/2) (71) 

n,m=-oo "'-"o 

where 

K = exp{z7r^((5mod l^/a} 



Uo{T) = exp{-z7rr/(Pmod ly/a}, (72) 
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so that U{T) = KUq. [ Here again we employ a convention where capital letters denote 
scaled variables i.e. P = p/b, Q = q/a, V = v/a and below X = x/a.] The matrix element 
in Eq. ( [7lD can be written as 



{q - v/2\kUoiT)fn,mU^oiT)k^q + v/2) = 

r dq' e^^'^^^'^q - v/2\kUo{T)\q' - naa/2){q' + naa/2\U^Q{T)k^q + v/2), (73) 

J —oo 

which involves matrix elements over kUQ{T). These matrix elements may be rewritten as 
(see Appendix B) 

oo 

(a;|^f7o(r)|x') = exp{i7r^(Xmod if /a} J2 9kiv)Six - ^' - kaa), (74) 

fc=— oo 



where 



dv e-^""'''/"e2^^*^^ (75) 

Jo 



9k[z, 

Inserting Eq. ([7^) into Eq. (|73|) gives 

(g - v/2\kUo{T)f^,rrM{T)k^\q + v/2) = 

oo 
fc,i=— oo 

g2.im(Q-(fc+/)a/2)^^^ - (n + / - k)aa). (76) 
Inserting Eq. (|7|) into Eq. (|7l|) gives, after some algebra 

-| oo 

p(p,g,r) = - Pn,™(o)e2-("^+'"«)«:„,„(p,g) (77) 

ao 



n,m=—QO 



with 

00 

Kn,m{P^Q)= E 9M9Me'^"'^'-'^'' ■ 
fc,/=— 00 

g-i7r«(fc+/)mgi7r5[(Q-a(n-fc+/)/2)mod l]'^ / a ^-iTT£,[{Q+a{n~k+l) /2)xiiod ^f^g'j 

Equations ( |77D and (|78D provide an expression for propagation of the quantum map in phase 
space. However, extensive algebraic manipulation is necessary to extract the propagator in 
the Fourier basis. We begin by noting that we can write 
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rb 

Pn,m 



dp J cigop(Po,Qo,0)e-2-("^o+'"«") (79) 

JO 

so that Eq. ( [77D can be written in the form 

p(P, Q, T) = t dpo r dqo p(Po, Qo, 0)/c(P, Q, T; Pq, Qo, 0) (80) 

JO JO 

where 

1 oo 

/s:(P,Q,T;Po,Qo,0) = - E e2-KP-Po)+-(Q-Qo)]^^^^(p^ g) (gi) 

n,m=— oo 

is the kernel of the Liouville propagator in the Wigner-Weyl representation. Since 

p{P, Q,T) = \y1 Pn™(r)e2-("^+-«) (82) 

n,m 

the Fourier coefficients of the propagated Wigner functions are 

Pn,miT) = f'dpo pdqo p(Po,Qo,T)e-2-*('^^«+"^««). (83) 

JO JO 

Given Eqs. (|80|) - (^21), the goal is to rewrite the propagation of Pn,m(0) in terms of Eq. 
(|65|), so as to extract G{n, m; k, I), i.e. in the first iteration we want 

oo 

Pn,miT)= E Gin,m;k,l)pkM- (84) 

fc,Z=— oo 

We begin by noting that Eq. (^Tj) can be manipulated, in accord with Appendix C, into the 



form 

1 



«:(P,Q,T;Po,Qo,0) = - E e^-I^^+^^^l ■ 

n,m=— oo 

oo 

k,l = — OD 

^iTTi[(Qo-a{n-k-l)/2)mod 1]"^ /a^-iTr^[(Qo+a(n+k+l)/2)mod Xfja ^gg-j 

and so Eq. (|80|) can be written as 

1 rb ra 

p{P,Q,T) = - E ^2.i[nP+mQ] ^^1 piPo,Q 0,0) 

ah „ J-'_^ Jo Jo 



n,m=—oo 



k,l=—oo 

pj7r?[{Qo-a(n-fc-0/2)mod l]'^/a -in^[{Qo+ain+k+l)/2)mod l]^/a 
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(86) 



Comparing Eq. (|^) to Eq. (|82|) we identify 

rb ra 

pn,m{T)= dpo dqo p{Po,Qo,0) 
Jo Jo 



k,l=—oo 

Jn^[iQo-oi{n-k-l)/2)mod l]'^/a -ini[{Qo+a{n+k+l)/2)mod l]^/a 



(87) 



Using the identity (see Appendix B for a proof) 

oo 

E ^7,(.)e-^-^ (88) 



oo 

-iTTz{xmod l)^/a \ ^ n ( ■y^ /d^^""*^-? 



j = -oo 

one obtains 

oo 

i,j,k,l=—OD 

p„+fc_z,^+,_,(0)e-^™('=+')(™+*-^)e-*™(*+^>. (89) 

This can immediately be generahzed to the iterative form giving pn,m{{M + 1)T) in terms of 
Pn,m at time MT. However, although this equation correctly describes the manner in which 
the Fourier coefficients of a Wigner distribution are mapped under the quantum Liouville 
propagator it is too complicated to be of use. Let k' = n + k — I and i' = m + i — j. Under 
this change in the integers of summation Eq. (^) becomes 

oo 

Pn,m{T) = Y 9i'+j-mi09*j{09k'+l~n{v)9!{v) " 
i' ,j,k' ,l=—oo 

Pk',i' (0)e~*™('='+2'-")*'e"'™(''+^^'"™)". (90) 
But the sums over j and I may be done explicitly. To begin with 

oo oo 1 1 

E ^/;(0^7.'+.-™(0e-''^^(°"^^' = Y J diyj du'. 

j=—oo j=~oo ^ ^ 

^iTT^i/'^ /a ^-2mju ^-in^u'^ /a ^2m(i'+j-m)iy' ^-2m{an)j ^q-^^ 

and use of identity ([T7|) yields 

oo „i „i 

j=-oo ^0 ^0 

g-i^5v'V"g2«(i'-m)^'^(^/ _ ^ an)mod 1) (92) 

^j^g— j7r5[(i^+an)mod 1]'^ /a^iTv^i/'^ /a^2TTi{i' —m){u+an) (93) 
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Similar calculations give 

l=-oo -^0 

Insertion of these results into Eq. (PO), setting k' = k and i' = I, and simplifying, gives the 



desired result, i.e. Eq. ( p4D with 

G{n, m- k, I) = e^™(^'-"'^) dv (;;y'e-^-«[('^+°«)>^°d i]V"ei-«^V« 



^■Ki(l—m)u ^—mri[(u' +al)mod 1] /a^i-KrjV /a^TTi{k—n)u' (95) 



as the propagator of the Fourier coefficients [p3]. This can be immediately generalized to 



give Eq. ( |65D with the same propagator G. The crucial difference between the classical Gc 
[Eqs.(^),(p5D] and the quantum propagator G is clear insofar as the quantum propagator 
mixes contributions from all k,l components to produce each n, m whereas Gc does not. 



Note that propagator (|95D moves distributions forward in time. To obtain its inverse, 
i.e. the propagator which moves distributions backward in time, simply take the complex 
conjugate of Eq. (|95| ) and interchange n with fc, and m with I (see Appendix D for a proof). 

As a simple check on the validity of the quantization procedure we redo the "free particle" 
case i.e. ^ = 0. In this case G is given by 



and the eigenequation is 

oo 

y G{n,m;k,l)pk,i = e-'^'^Pn,m- (97) 



A;,Z=— oo 

-ike, 



Consider a solution of the form p^ i = fi^e Substitution and use of Eq. ([T7|) yields 



-inamn 1^ ^^^-i^^[(v+am)mod \f j ^ ^mv^ j ^ ^-2i.inv _ _ am/2)mod 1) = e-^^^C"^"^" 

27r 



(98) 



and this implies that 



A = ^{[(— + «m/2)mod 1]^ - [(^ - «m/2)mod 1]^}. (99) 
aT 27r 27r 

Choosing ^ = ^ where po G [0, b) gives the same solutions [Eqs. (^) and (^)] which we 
obtained previously. 
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V. CORRESPONDENCE 



Given Eq. (^) it is straightforward to examine the dynamics in the classical limit 
(a — i> 0). We do so at fixed finite time so that, with respect to the issue of ergodic properties, 
we are taking the correct order of limits, /i — > prior to any long time limit 0. 

For each set of values {u, u', n, 1} there exists a sufficiently small a such that [u + 
an)mod 1 = u + an and (z/' + a/)mod 1 = u' + al. Thus in the a — >• limit 

Jo Jo 

= S{l,k),^'in,m) (100) 

which is the correct classical limit, Eq. (^). Note that this result holds for all maps of the 
form given by Eq. (0), including the classical chaotic cat map. 

Having demonstrated that these quantum maps give the correct classical limit it is nec- 
essary to obtain its dynamical characteristics in the case of h ^ 0. It is straightforward to 
demonstrate that, independent of the value of r], ^, the quantum propagator has at least two 
eigenf unctions with unit eigenvalue and hence that the quantum dynamics is nonergodic. 
Two such eigenfunctions are the uniform distribution and the Schrodinger propagator. To 
see this consider the uniform distribution lfc,z(0) = SkfiSi^ under propagation. Using Eq. 
(p5|) we have 

(T) = G{n,m;0,0) 

^—2TTimna I j ^—iTT^[{u+an)mod 1]'^ /a ^itt^l''^ /a ^—2TTimi/ ^—2TTinu' 

Jo Jo 

= (5„.o f'du 
' Jo 

= Snfi^mfi- (101) 

Thus the uniform distribution is an eigenfunction with unit eigenvalue. [This result also 
implies that the dynamics is measure (i.e. area) preserving.] 
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Second, note that the Schrodinger propagator is an eigenfunction of the Liouville propa- 
gator since tj'^{U)tj = tj^^{U)JJ = U. It remains to demonstrate that the Wigner function 
associated with U is on a unit celL To see this note the form of U in the Wigner- Weyl 
representation: 

CO 

U('"){P,Q)= J2 ^-n(r7)c/;;(Oe"'""""e2"*("^+™Q\ (102) 



n,m=— CO 



Then 



1 rl rl rl <^ 

dP dQ\U^^\P,Q)\'= dP dQ 9-kiv)9ii09-n{v)9*mi0 " 

Jo Jo Jo 77 

n,m,k,l=—oo 

^2ni[{n—k)P+{m—l)Q]^—iTT{nm—kl)a (103) 
OO oo 

= E \9k{v)\'- E \9im' (104) 

k=—oo l=—oo 

= 1^ = 1 (105) 



smce 



E 9j+k{z)9*k{z) = 6j,o- (106) 

A;=— oo 

The latter identity is proven in Appendix B. 

Thus is an eigenfunction of the mapping with eigenvalue 1, as is the uniform 



distribution and so |]18[ the mappings are not ergodic. Note, furthermore, the a — > limit 
of U^'^'^ displays an essential singularity. Hence this eigenfunction does not exist classically, 
a result which is essential to achieve the proper ergodic classical limit for classically ergodic 
cases like the cat map. 

The extended quantization procedure therefore provides a straightforward treatment of 
the cat map, a system which has been the subject of considerable controversy. 



VI. WHY NO HAMILTONIAN EIGENFUNCTIONS? 



In addition to the explicit results on maps, this work has also introduced a quantiza- 
tion procedure leading to results which differ significantly from the traditional Schrodinger 
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equation approach. This difference arises primarily from the fact that periodic boundary 
conditions apphed to the density matrix [T^^pT = p for T = Tq{a) and T = Tp(6)] allow 
for solutions when similar boundary conditions, applied to the wavefunction, preclude the 
possibility of a solution. Formally this is so because periodic boundary conditions on the 
wavefunctions (corresponding to the condition Tp = p, pT = p, on the density matrix) is a 
far more restrictive requirement than the condition T~^pT = p. Nonetheless it is of interest 
to examine the character of the mapping dynamics and to expose the physics underlying the 
lack of wavefunctions. 



Consider the Arnold cat map [T^ as a specific example. Written explicitly, the cat map 
(^7 = ^ = 1 ill Eq. (0) reads 

Q„+imod 1 — Qninod 1 = P„mod 1 + V{Qnin^od 1, P„mod 1) (107) 

P„+imod 1 — P„mod 1 = P^mod 1 + Qninod 1 + P((5„mod 1, P„mod 1) (108) 

where 

V{Qmod 1, Pmod 1) = -IxA^^uA^iQraod 1, Pmod 1), (109) 

P(Qmod 1, Pmod 1) = -IxAiUA^iQ^od 1, Pmod 1) - 2xA4(Qmod 1, Pmod 1). (110) 

Here xb is the characteristic function on the set B and the phase space regions Ai are defined 
by: 

A2 = {^< Qmod 1 < 1} X {1 - 2gmod 1 < Pmod 1 < 1 - Qmod 1} (111) 
As = {0 < gmod 1 < i} X {1 - Qmod 1 < Pmod 1 < 2 - 2gmod 1} (112) 
^4 = < Qmod 1 < 1} X {2 - 2Qmod 1 < Pmod 1 < 1}. (113) 



Equations ( |107| ) and ( |108|) show that the changes in q and p under the mapping depend upon 



an external boost V and external force F which are functions of position and momentum. 
Thus the cat map equations are a discretization of a Langevin type process 

— =p/p + V{q,p) 

|._cp-|£ + ffe„ (114) 
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with a negative coefficient of friction ( = —1 and potential $ = — Since this set of 
equations describes the Hamiltonian dynamics of a particle moving in the vicinity of a single 
hyperbolic point ifV = ( = F = 0, Eqs. (44) and (45) constitute a set of perturbed dynamics 
around an unstable point. Thus the classical cat map can be regarded as a system subject to 
external boosts and forces, rather than a stroboscopic viewing of a true Hamiltonian system. 
Since the system is one which is subjected to external boosts and forces, it does not admit 
a pure state description. 

VII. SUMMARY 

In summary, we have introduced a von Neumann based quantization procedure in concert 
with Bloch type boundary conditions in order to achieve a physically complete quantization 
of mappings on a torus fl^''. The more traditional approach [^,0 of quantizing the mappings 
via the Schrodinger equation with associated periodic boundary conditions yields solutions 
only for tori satisfying ab = hN, and excessively restricts the mapping dynamics to the 
dynamics on a set of rational points in {q,p) space, associated with a discrete Liouville 
spectrum. Alternative approaches such as applying Bloch type boundary conditions on 
the wavefunction or periodic boundary conditions on the density matrix encounter similar 
difficulties. From the perspective of the density matrix periodic boundary conditions imply 
Tp = p,pT = p for T = Tq{a) and T = Tp{b). Adopting Bloch boundary conditions on 
the density matrix {T^^pT = p), the only remaining possibility, allows for a broader class 
of solutions obtainable even when wavefunctions do not exist. The fact that pure states do 
not in general exist is a consequence of the fact that boundary conditions here correspond 
to the application of external forces to the system. 

The resultant quantum system has a continuous Liouville spectrum, and the von Neu- 
mann quantized cat map dynamics is not ergodic until the proper classical limit is reached. 
This result provides the ffist demonstration of a straightforward limiting procedure which 
completely recovers the full chaotic classical dynamics in the h ^ limit. Further studies 
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are in progress to ascertain the character of the quantum dynamics and its dependence on h 
as the system approaches the classical limit. Further, this study suggests the possibility of 
quantization via a generalized von Neumann (or continuity) equation for systems described 
as a Langevin type process i.e. Eqs. ( |114| ), a fertile general area for further study. 

Acknowledgment: Support from the Centre of Excellence in Molecular and Interfacial 
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APPENDIX A 



Here we discuss briefly the algebra generated by the operators fn,m- It is triv- 
ial to show that these operators satisfy the following conditions; (1) closure fn,mfk,i = 
e*''^"'""'')°/„+fc,m+i, (2) associativity {fn,mfk,i)fi,j = fn,m{fk,ifi,j), (3) the existence of an 
identity /o,o = 1, (4) the existence of inverses fn,mf-n-m = 1 = f~n-mfn,m, and (5) 
fnm = f-n-m = {fn,m)~^ which implies that they are unitary operators. Furthermore, 6) 
every element of this algebra can be generated from the set {/i,o, f-i,o, fo,i, /o,-i}- 



APPENDIX B 

Here we will prove some of the simple formulae used in this paper. Underlying all of these 
formulae is the identity in Eq. ([T7|) . We have defined integrals gk{z) = dve~'^'^^^^ /°'e^'^'^^^ 
and here we will prove two simple identities and one formula specific to the quantum map 
problem involving these integrals. 

The first identity is that of Eq. (^) . Using identity 



Jo 

„l oo 



^■Ki{v—x)j^—i-KZv'^/a 



/•i ^ 

j=-oo 
oo 

= gA^y-^^'"''- (115) 

j=-oo 

The second identity is that of Eq. (|106|) . Again using identity ([T7| ) 

oo ^\ oo 

• - ^0 ^0 j^oo 

Jo 

The formula specific to the quantum map problem is Eq. (0). The proof is as follows; 
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{x\kUo(T)\x') = / dp{x\p){p\x')e-'^''^^'^°^ ^^'/^ 



oo 



gi7r5{Xmod if /a I ^ ^ipjx-x') /h ^-iiTr](Pmod 1)2/q 



oo 



27Tn 



— dp 



^i7r5(Xmod l)^/a ST^ ( \ / "/^ ip{x-x'-aak)/h 



k=—oo 
oo 

fc=— oo 



6 ' 
oo 277/1 



;=— 00 
00 

gk{r])6{x - x' - aak), (117) 

;=— 00 

where we have used identity (^) and the closure relation for momentum. 

APPENDIX C 

Here we will show that Eq. (^) can be put in the form of Eq. (pS]). To begin, we 
substitute expression (|7|) for Kn,m{P,Q) into Eq. (|8TD. This gives 

00 

^2ni[n{P-Po)+m{Q-Qo)] 



n,m=— 00 



J2 9k{,r])g*i{vi)e^ 

k,l = — OD 

'iTTa(k+l)m in^[{Q-a{n-k+l)/2)mod 1]^ /a -in^[{Q+a{n~k+l) /2)mod l]^/a 



KiP,Q,T;Po,Qo,0)= E e 

n,m 

^*^^ymii-k)p 

^~tTTa(k+l)m ^nr^[{Q-a{n-k+l)/2)mod /a ^-nT^[(Q+a(n~k+l) /-Zjaiod l\^/a (118) 

Using identity (^^ for the sum over m then gives 

oo 

«:(P, Q, T; Po, Qo, 0) = E e2™(^-^")5(g - Qo - (A; + l)a/2 - m) ■ 

•n,m,k,l=—OD 

g^^^-Sjg*^^^^27Tiil-k)P^in(l{Q-ain-k+l)/2)mod l]"" /a^-i7Ti[iQ+a(n-k+l)/2)mod ifja ^^-yO)) 
oo 

= E e2™(^-^«)5(g - - (A; + /)a/2 - m) ■ 

n,m,fc,i=— oo 

^^(•^)^*(^^)g27rj(/-A:)Pgi7r?[(Qo-a(n-2A:)/2)mod l]Vag-i7rC[{Qo+a(n+20/2)mod (^]^20) 

Now using identity ([TtD in reverse we obtain 

oo 

«:(P, Q, T; Po, Qo, 0) = E e2"["(^-^°)+"^(«-«°)l(?,(r/)(?r(r^)e2^^('-^)^ ■ 

n,m,fc,i=— oo 

g-i7ra{fc+Omgi7r^[(Qo-Q(«-2fc)/2)mod l]2/Q^-j7rg[{Qo+a{n+2/)/2)mod (121) 

Making the change in the integer of summation n' = n + 1 — k we then obtain 
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H{P, Q, T; Po, Qo, 0) = ^ ^2ni(n'P+mQ)g^^^^g*^^•^^-2.^[(n'+k-l)Po+mQo] . 

n' ,m,,k,l=—oo 

A:-0/2)mod 1]^ /a ^-i7T^[{Qo+a{n' +k+l) /2)mod (122) 

Now setting n' = n and changing the order of the factors we obtain Eq. (^). 

APPENDIX D 

Here we will show that taking the complex conjugate of propagator ( pSj) and interchang- 
ing n and fc, and m and / gives the inverse of propagator (^). This inverse propagator 
propagates distributions backward in time. The inverse propagator according to this pre- 
scription is 

Jo Jo 

^2iTi{l—m)i' ^iTTTiKu' +am)mod l]'^ /a ^—inriu''^ /a ^2ni(k—n)u' (123) 

To verify that is indeed the inverse of G we need only show that 

oo oo 

G{n,m]r,s)G'^{r,s;k,l) = 6n,kSm,i = G'^{n,m;r, s)G{r, s; k,l). (124) 

r,s=— oo r,s=—OD 

We begin with the first equality; 

oo CO 1 1 I I 

y G{n,m;r,s)G-\r,s;k,l) = e'''^^^-''"'^'' V dv dv' dv" dv'" ■ 
„_ „ „„_ „Jo Jo Jo Jo 



r,s=—oo 



—iTT£^[{u+an)mod iTr^i/^/a 27ri(s— — i7rr;[(f'+os)mod 1]'^ /a inrju''^ / a 2TTi(r—n)u' 

in^Kiy" +ak)mod 1]"^ /a -in^u''^ /a 27Ti{l-s)i^'' i7vr][{u''' -\'as)mod 1]"^ /oi -inrji'"''^ /oi 2'Ki{k-r)v"' /-i r)r\ 

c o o o o c . I 

Using identity ([T7| ) to do the sum over r gives the delta function 5{v"' — v') and then doing 
the integral over v'" we obtain 

oo /.I fl 

V G(n,m;r,s)G-^(r,s; = e^"^'^'-"™)" V dv dv' dv" ■ 

r,s=-oo r,s=-oo-^0 ^0 ^0 

— i7r5[(i^+an)mod l]^/a inS^u^/a 2TTi{s—m)i' —iTTr^[{u'+as)iaod l]'^/a i-K'qu'^/a 2TTi(r-—n)u' 

^i-K^[{u" +ak)mod 1]'^ /a ^—iiT^i'"'^ /a ^2ni{l—s)u" ^inriKu' +as)mod 1]'^ /a ^—i-wqu'^ /a^2ni{k—r)i>' (126) 
oo „i „i „i 

= e*"('='-""^)° V dv dv' dv" ■ 

r,s=-oc-JO Jo Jo 

— 27r^[(i^+Q;n)mod 1]'^ / a Atv^i^'^ / a 2ni{s—rn)u 2ni{k~n)iy' i7v^[(i/'' -\-ak)mod 1]'^ /a —iir^u"^ /a 27vi{l — s)v'" 07^ 
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Now performing the sum over s, again through the use of identity (p!?!), gives the delta 
function (5(z/" — u) and doing the integral over u" we obtain 



V G{n, m; r, s)G-\r, s; k, I) = e^'^^'^'""™)" / du / (i;y'e-*-«[(^+"n)mod . 
r,«=-oo -^0 -^0 



Doing the integral over gives a factor Sn,k and we obtain 

„_ „ Jo 



r,s=— OO 



^27ri(Z— m)j^gi7r5[(j/+afc)mod i7rgi/^/a^ ^ (129) 



= e'<kl-nm)a dv e^'^'^^''^^ 5 n k ■ (130) 
^0 

Finally, doing the integral over v we obtain a factor 5m,i and so 

^ G{n,m;r,s)G~^{r,s;k,l) = Sn,kSm,i (131) 

r,s=— OO 

as we claimed. The identity in Eq. (|124|) can be proven in the same fashion. 
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